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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2145 


LIFT- CANCELLATION  TECHNIQUE  IN  LINEARIZED 
SUPERSONIC-WING  THEORY 
By  Harold  Mire Is 


SUMMARY 

A  lift- cancellation  technique  is  presented  for  determining 
load  distributions  on  thin  wings  at  supersonic  speeds.  The  tech¬ 
nique  retains  certain  features  of  the  method  recently  introduced 
by  Theodore  R.  Goodman,  while  simplifying  and  generalizing  others. 

A  general  expression  is  derived  for  the  load  distribution  over 
a  cancellation  wing.  This  expression  permits  the  determination  of 
lift  distributions  on  wings  that  cannot  be  solved  by  cancellation 
techniques  based  on  the  superposition  of  conical  flows.  The  bound¬ 
ary  conditions  for  either  a  subsonic  leading  edge  or  a  subsonic 
trailing  edge  can  be  satisfied.  Applications  of  the  expression  to 
swept  wings  having  curvilinear  plan  forms  and  to  wings  having 
reentrant  side  edges  are  indicated. 


INTRODUCTION 

The  method  of  lift  cancellation  for  obtaining  the  lift  dis¬ 
tribution  on  thin  wings  at  supersonic  speeds  was  first  suggested 
in  reference  1.  The  lift  distribution  on  a  given  wing  is  deter¬ 
mined  by  canceling  excess  lift,  through  the  use  of  a  "cancellation 
wing,"  on  a  related  plan  form  having  a  known  loading.  This  approach 
has  been  applied  by  several  authors  (for  example,  references  2  to  4). 
The  expressions  provided  in  reference  1  are  applicable  for  wings 
that  can  be  generated  by  the  superposition  of  conical  fields. 

A  procedure  is  presented  in  reference  5  for  determining  lift 
on  a  more  general  class  of  plan  forms  than  can  be  handled  by  coni¬ 
cal  superposition.  The  method  utilizes  a  surface  distribution  of 
doublets  and  an  inversion  by  means  of  Abel's  integral  equation  and 
is  equivalent  to  a  lift  cancellation. 
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This  report ,  prepared  at  the  NACA  Lewis  laboratory,  retains 
certain  features  of  reference  5  (that  is,  the  use  of  a  surface  dis¬ 
tribution  of  doublets  and  an  inversion  by  means  of  Abel's  integral 
equation),  whereas  other  features  are  simplified  and  generalized. 

The  simplification  consists  in  eliminating  steps  in  the  procedure 
for  obtaining  lift  distributions.  The  generalization  consists  in 
determining  a  solution  that  can  be  made  to  satisfy  the  boundary 
conditions  for  either  a  subsonic  leading  edge  or  a  subsonic  trailing 
edge  (Kutta  condition)  •  The  method  of  reference  5  yields  only  the 
Kutta  solution.  The  lift -cancellation  technique  developed  herein 
is  illustrated  by  several  examples . 

In  a  concurrent  investigation  (reference  6) ,  source  distribu¬ 
tions  and  integral-equation  formulations  have  been  applied  to  obtain 
the  loading  on  a  special  series  of  cancellation  wings.  Reference  7 
employs  some  of  these  cancellation  wings  for  the  determination  of 
lift  and  moments  on  swept  wings. 


THEORY 

The  usual  assumptions  of  an  inviscid  fluid  and  small  perturba¬ 
tions  are  made.  The  velocity  field  consists  of  the  free-stream 
velocity  U  (taken  in  the  positive  x-direction)  plus  the  perturba¬ 
tion  velocities  u,  v,  and  w.  The  wing  boundary  conditions  are 
specified  in  the  z  =  0  plane. 

The  local  lift  coefficient  ACp  may  be  expressed  in  terms 
of  Au.  That  is. 


_  PB-PT  _  2(uT"ub)  _  2Au  f-n 
P  q  U  _  U 

(All  symbols  used  in  this  report  are  defined  in  appendix  A.)  Inas¬ 
much  as  the  local  lift  coefficient  is  directly  proportional  to  Au, 
Au  will  be  referred  to  as  "lift"  in  later  developments. 


Lift-Cancellation  Method 

The  lift  distribution  on  a  given  wing  is  to  be  determined  by 
canceling  excess  lift  on  a  related  wing  with  a  known  loading.  The 
method  is  illustrated  in  figure  1.  The  wing  for  which  the  lift 
distribution  is  desired  is  shown  in  figure  1(a)  The  solution  can 
be  expressed  as  the  two-dimensional  wing  (fig.  1(b))  minus  a  can¬ 
cellation  wing  (fig.  1(c)).  The  loading  in  region  I  of  the 
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cancellation  wing  equals  the  loading  in  the  corresponding  region  of 
the  two-dimensional  wing  and  the  upwash  w  in  region  II  of  the 
cancellation  wing  is  zero.  The  loading  of  the  two-dimensional  wing 
minus  that  of  the  cancellation  wing  satisfies  the  boundary  conditions 
for  the  flow  about  the  given  wing  and  is  the  desired  solution. 

The  fundamental  problem  in  the  lift-cancellation  method  is  then 
to  determine  the  lift  in  region  II  of  a  cancellation  wing  subject 
to  the  condition  w  =  0  in  this  region  and  with  the  assumption  of 
a  known  loading  in  region  I.  Solution  of  this  problem  is  presented 
in  the  following  sections. 


Derivation  of  Lift-Cancellation  Equations 

The  lift  distribution  in  region  II  will  be  expressed  in  terms 
of  quantities  in  region  I. 

Consider  the  cancellation  wing  shown  in  figure  2.  The  portion 
of  the  leading  edge  to  the  left  of  the  origin  coincides  with  a  Mach 
line.  The  portion  of  the  leading  edge  to  the  right 
(designated  r  =  r^(s))  is  shown  as  a  supersonic  edge,  although  no 
restrictions  as  to  a  subsonic  or  supersonic  edge  are  imposed.  (A 
plan- form  edge  is  subsonic  or  supersonic  depending  on  whether  the 
component  of  the  free  stream  normal  to  the  edge  is  subsonic  or 
supersonic.)  The  line  designated  r  =  ^(s)  separates  region  I 
and  region  II  and  is  assumed  to  be  subsonically  inclined  to  the 
free  stream  at  all  points.  This  line  corresponds  to  a  plan-form 
edge  of  the  wing  for  which  the  lift  distribution  is  desired. 

General  solution  for  load  distribution  on  cancellation  wing.  - 
The  upwash  field  in^the  z  =  0  plane  (due  to  an  arbitrary  distribu- 
tion  of  vorticity  A u  and  Av)  may  be  written,  from  reference  8, 
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The  symbol  |  designates  the  finite  part  of  an  infinite  integral, 
as  defined  in  reference  9.  Application  of  the  finite-part  concept 
to  linearized  supersonic-wing  theory  and  the  evaluation  of  the  finite 
part  of  an  infinite  integral  are  discussed  in  references  8  and  10. 

For  the  present,  it  will  suffice  to  state  the  fundamental  definition 
of  the  finite  part  of  an  integral  with  a  3/2-power  singularity, 
namely. 


By  a  transformation  to  the  Mach  coordinates  of  reference  11, 


x  =  I  (s+r)  r  =  (x-Py) 

M  op 

y  =  |(s-r)  s=^(x+py) 


elemental  area  = 


dr  ds 


equation  (2)  becomes 


Upon  substitution  of  the  limits  of  integration,  as  indicated  in 
figure  2, 
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Integrating  by  parts,  noting  that  A<p  =  0  at  rQ  =  r^(s0),  and 
recalling  the  definition  of  the  finite  part  (equation  (3))  yield 


(9) 
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The  right  side  of  equation  (14)  vill  he  considered  known.  Equa¬ 
tion  (14)  is  then  an  integral  equation  for  A<P-q.  The  solution 
(appendix  B)  is 


i  2V  o' 
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Equation  (15)  indicates  that  the  doublet  strength  in  region  II, 
namely  A<Pjj,  can  obtained  by  a  line  integration  along  sQ  «  s 

in  region  I,  The  geometric  interpretation  of  the  various  terms  in 
equation  (15)  is  shown  in  figure  3. 

It  can  be  shown,  by  expanding  $Pj  in  a  Taylor’s  series  about 
r  =  r2(s),  that  equation  (15)  yields  a  continuous  solution 
(ACpjj  =  A^pj)  at  r  =  r2(s).  (A  discontinuity  in  A<P  implies  a 
lifting  line  (reference  12)  and  is  unrealistic.) 


The  lift  distribution  in  region  II  can  be  expressed  as 


Auii  “ 


dx 


5AtpII  Br  ^AqpII  Bs 

Br  SEE  Bs  chE 


or,  from  equation  (15), 
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Differentiation  yields  (See  appendix  C.) 

_ rr2(s) 


Au 


r-r2(s) 


II 


AuT  dr 

1 _ 2. 


K)#K 


1  - 


dr2(s) 


nr9(s) 


ds 


2P«/^r-r2(s) 


(8AuI-AvI)dro 


(s) 


A/r2(s>" 


(17a) 


Equation  (17a)  is  the  desired  expression  for  the  lift  distribution 
in  region  II  in  terms  of  quantities  in  region  I. 

Consider  Au^j  to  consist  of  two  components,  Aujj*  and 
Au^",  where  Au^'  and  Au^,f  are  the  and  second  terms 

on  the  right  side  of  equation  (17a),  respectively.  Investigation 
of  the  integrals  indicates  that  at  r  =  r2(s),  Aujj*  =  Auj; 

whereas  Aujj",  in  general,  has  a  half-order  singularity. 

When  region  II  is  to  the  right  of  region  I  (fig*  3(b)),  the 
integration  for  Aujj  is  conducted  along  the  line  rQ  =  r  and 

may  be  written  as 

^s2(r) 
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Discussion  of  equations  (17a)  and  (17b).  -  In  the  paragraph 
preceding  equation  (2),  the  line  r  =  r2(s)  was  described  as  sub- 

sonically  inclined  at  all  points  to  the  free  stream.  This  condition 
is  necessary  so  that  the  inner  integral  in  equation  (lla)  (that  is, 
G(r,sc))  can  be  equated  to  zero  for  all  points  in  region  II.  If 

this  restriction  on  r  =  r2(s)  is  not  satisfied,  the  development 
beyond  equation  (lla)  becomes  invalid.  The  derivation  of  cancella¬ 
tion  equations  when  r  =  r2(s)  is  supersonic  was  not  undertaken 

because  such  problems  can  be  solved  more  simply  by  other  methods. 

In  regard  to  the  boundary  conditions,  it  has  been  assumed  that 
AUj  is  specified.  Equations  (17a)  and  (17b),  however,  indicate 

that  a  knowledge  of  Avj  is  also  required  in  order  to  obtain  a 
solution  for  Aujj.  Two  possibilities  exist,  as  illustrated  in 

figure  4.  In  the  first  case  (fig.  4(a)),  region  I  is  upstream  of 
region  II  (along  the  line  r  =  r2(s))  and  Avj  is  uniquely  defined 

by  the  specified  Auj  according  to  the  relation 

AvI  =  j  AuI  <^0  (18a) 

yJxL(y) 

The  integration  is  conducted  along  lines  of  constant  y.  In  the 
second  case  (fig.  4(b)),  region  II  is  upstream  of  region  I  (along 
the  line  r  =  r2(s))  and  the  expression  for  Avj  in  region  I 

(for  y  <  0)  is 


Equation  (18b)  indicates  that  a  knowledge  of  Aujj  is  required  in 
order  to  find  Avj.  But  Avj  must  be  known  (equation  (17a))  before 
Aujj  can  be  found.  Thus,  the  solution  for  Aujj  from  a  specified 
Auj  is  not  unique  for  the  configuration  of  figure  4(b)  and  an 
additional  boundary  condition  must  be  imposed.  The  line  r  =  r2(s), 
however,  corresponds  to  a  plan-form  edge  of  the  airfoil  whose  load 
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distribution  is  desired.  The  situation  indicated  in  figure  4(b) 
occurs  vhen  r  =  r2(s)  corresponds  to  a  subsonic  trailing  edge. 

The  additional  condition  to  be  imposed  is  therefore  the  Kutta  con¬ 
dition.  In  terms  of  the  cancellation  wing,  this  condition  requires 
that  the  perturbation  velocities  be  continuous  in  crossing  r  =  r2(s). 

Solution  for  Aujj  satisfying  Kutta  condition  at  r  =  ^(s).  - 
It  will  now  be  shown  that  when  the  Kutta  condition  is  imposed  at 
r  =  r2(s),  the  appropriate  Avj  distribution  is  such  as  to  make 

the  second  integral  in  equation  (17a)  identically  zero;  that  is. 


pr2(s) 


^(s) 


gAUj-Av- 


[r2(s)-ro] 


172  =  ° 


dA<T>i 

or,  inasmuch  as  0Auj-Avj  =  M  — — , 


pr2(s) 


'SrT  dr< 


i-l(s) 


[r2(s)'ro]l/2 


=  0 


This  concept  and  its  proof  follow  from  a  suggestion  of  H.  S.  Eibner 
of  the  NACA  Lewis  laboratory. 


Thuf,  from  equation  (7),  (12),  and  (13), 
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=  0 


for  all  points  (r,s)  in  region  II.  Therefore, 


,r~(s) 


rl(s) 


dACpx 


(r-rj 


lj2 


r>r 


r2(s) 


(r-rQ) 


172 


(19) 


(20) 


t 


« 


1297 


MCA  TN  2145 


11 


Upon  taking  the  limit  as  r  approaches  r2(s) ,  equation  (20) 
becomes 


^r2(s) 

rX 

Ba<Pt 

SACPjj 

Sl*°,  -  lln 

ri(s)  -**(•> 

0 

r2(s)  J 

(21) 


However,  <5t$Pj/br0  must  be  continuous  in  the  vicinity  of  r2(s). 
(The  perturbation  velocities  on  the  basic  wing  can  be  discontinuous 
only  along  Mach  lines  or  along  plan-form  edges.  Inasmuch  as 
r  =  r2(s)  is  neither  of  these  cases,  all  derivatives  of  ACpj  must 
be  continuous  in  the  vicinity  of  r  =  r2(s).)  When  the  Kutta  con¬ 
dition  is  imposed,  dACpjj/dr0  is  therefore  also  continuous  (and 
bounded)  in  the  neighborhood  of  r  =  r2(s).  Then,  using  a  mean 
value  for  dACp jj/dr0, 

rX 
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=  0 


Therefore 


(22) 


pr2(s) 


rl  (s) 


^7  6X0 


[r2(s)-r0]1/2 
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(23) 


which  was  to  be  proved* 
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The  solution  for  Aujj  that  satisfies  the  Kutta  condition  at 
r  =  r2(s)  is  then,  from  equations  (17a)  and  (23), 


,r2(s) 


(s) 


Aui  drQ 

(r-ro)J^FT0 


(24a) 


for  the  wing  of  figure  3(a).  Similarly, 


s-sc>(r) 


>So(r) 


Mr) 


AuT  ds 
I  o 


(s-s0)A/s2(r)-s0 


(24b) 


for  the  wing  of  figure  3(b). 

An  alternate  derivation  of  equations  (24a)  and  (24b)  (appen¬ 
dix  D)  indicates  that  only  solutions  satisfying  the  Kutta  condition 
will  result  from  the  integral  equation  formulations  of  reference  5. 

Sidewash  in  region  II.  -  An  expression  for  Avjj  can  be  obtained 
by  differentiating  equation  (15)  with  respect  to  y.  The  result  is 


dr2(s) 
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Similarly,  for  region  II  to  the  right  of  region  I, 
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When  the  Kutta  condition  applies,  these  equations  become, 
respectively. 


It  should  be  noted  that  when  r  =  r2(s)  corresponds  to  a  subsonic 
trailing  edge,  Avj,  as  well  as  Avjj,  is  not  generally  known. 

The  preceding  expressions  are  therefore  primarily  useful  for  those 
problems  where  r  =  r2(s)  corresponds  to  a  subsonic  leading  edge. 


APPLICATIONS 

The  loading  in  region  II  of  a  cancellation  wing  is  given  by 
the  line  integrals  of  equation  (17a)  or  (17b),  When  the  Kutta 
condition  is  imposed  at  a  subsonic  trailing  edge,  the  expressions 
reduce  to  equations  (24a)  and  (24b).  These  equations  can  be  used 
to  find  the  load  distribution  on  a  large  variety  of  wings  that 
cannot  be  solved  by  cancellation  techniques  based  on  conical 
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superposition.  Wings  vitn  curvilinear  plan  forms  or  arbitrary 
camber  are  examples.  In  each  case,  however,  the  solution  for  the 
related  wing  must  be  known. 

The  equations  are  applied  in  several  illustrative  examples. 
Only  the  solution  associated  with  the  cancellation  wing  is  con¬ 
sidered.  The  complete  solution  consists  in  the  loading  of  the 
related  wing  minus  the  loading  of  the  cancellation  wing. 


Leading-Edge  and  Side-Edge  Cancellations 

In  these  cases,  the  lift  to  be  canceled  is  upstream  or  to  the 
side  of  the  plan  form  for  which  the  loading  is  desired  (figs.  5 
and  6). 

Tip  region  of  swept  wing.  -  The  loading  in  the  region  influ- 
enced  by  the  side  edge  (lla  and  IT^  of  fig.  5)  of  a  swept  wing 
having  a  subsonic  leading  and  a  supersonic  trailing  edge  can  be 
obtained  by  canceling  excess  lift  on  a  triangular  wing.  The  Kutta 
condition  is  applied  across  the  portion  of  r  =  r2(s)  influencing 

region  11^,.  The  lift  to  be  canceled  in  region  I  is  (reference  13, 
equation  (23)) 


AuI  = 


H02x 


H02(s+r) 


AJe2x2-fi2y 2  /^62(s+r)2-(8-r)2 


(25) 


where  H  and  9  are  constants  defined  in  appendix  A.  The  doublet 
distribution  in  region  Ia,  again  from  reference  13,  is 


A<P  =  H  ^02x2- p2y2  «  £5  y^02(s+r)2-(s-r)2 
Ia  M 

from  which 

AvT  .  ,  ___^_Ms-r) .  (26) 

a  y  A/e2x2-p2y2  A/e2(s+r)2_(s_r)2 

The  sidewash  distribution  in  region  I-jj  (that  is,  Avj^)  could  be 

found  by  an  integration  of  the  type  indicated  in  equation  (18b).  A 
knowledge  of  Avj^,  however,  is  unnecessary  in  the  present  problem 

because  the  Kutta  condition  is  applied  for  region  11^. 
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The  loading  in  region  IIa  is  obtained  by  substituting  equa¬ 
tions  (25)  and  (26),  with  r  replaced  by  rQ,  into  equation  (17a), 
which  yields 


Aun 
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H02  A^r-rjjTs") 

jt 


pr2(s) 

(s+r  )  dr 

_ v  o/  o _ 

r  (g)  (r-ro)^r2(s)_ro'\/e2(s+ro)2‘(s‘ro)2 
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02(s+ro)+(s-ro)]dr< 
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A/r2(s)_ro  AJe2(s+r0)2-(s-r0)2 


(27) 


For  region  II-^,  the  Kutta  condition  applies  and 


H02aJ 

pr2(s) 

(s+ro) 

1  drQ 

fr-r2(s) 

1 

K 

rl(s)  (r-r°H 

^2(s)-r0  AJ 

le2(s+ro)2-(s-r0)2 

(28) 

Equations  (27)  and  (28)  reduce  to  elliptic  integrals  of  the  first, 
second,  and  third  kind  upon  transforming  the  variable  of  integra¬ 
tion  from  rQ  to  C0Q  according  to  the  relation 

ro  =  I7e  Q1"0)®  +  a2®o2]  (29) 

where 

a2  =  (l+e)r2(s)-(l-e)s 
Equations  (27)  and  (28)  may  then  be  written 
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where 


a  =  (1+6)  r-(l-0)s 


a1  =  (1+0)  r1(s)-(l-0)s 


(31) 


St 

n  = - -  =  (1+0)  i*2(s)-(l-0)s 

Beentrant  side  edge.  -  A  plan  form  has  a  reentrant  edge  if  a 
line  of  constant  y  intersects  the  plan  form  at  more  than  two 
points. 

The  load  distribution  in  the  region  influenced  by  the  reentrant 
side  edge  is  to  be  determined  for  the  wing  of  figure  6.  The  side 
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edge  is  first,  for  simplicity,  the  straight  line  r  =  K2s,  which 
is  a  subsonic  trailing  edge  across  which  the  Kutta  condition  is 
applied.  The  side  edge  then  alternately  becomes  a  subsonic  leading 
and  a  subsonic  trailing  edge.  The  load  distribution  in  region  I  is 

simply  the  Ackeret  value  Au-j-  =  ^2*1,  and  Avj  =  0.  Regions  IIa, 

P  & 

lib,  and  II c  are  considered  separately. 

Region  IIa. 

From  equation  (24a)  with  AuT  =  ^2^ 

p 

K2s 

2aU  drQ 

s  e(r-ro)^K2s-r0 


Au 


II. 


.  H 


r-K2s 


=  MI  tan-1 

3« 


(K2+1)b 

r-K2s 


(32a) 


or,  in  x,y-coordinates 


4dU 


x+gy 


AUlIa  “  3*  tan  ^  p(m2x-y) 


(32b) 


Region  11^. 


A  knowledge  of  Av-^  is  required.  From  equation  (18b), 


=  3y 


4dU 

3* 


>y/m2 


ox 


tan' 


3y 


dx. 


y/®2 


2aU 


=  2aU 


V  k2-i 


(33) 
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A  cancellation  wing  may  induce  lift  that  itself  must  he  canceled  in  order  to  satisfy 
boundary  conditions  completely.  Thus,  in  figure  8,  the  cancellation  of  lift  in  region  I 


22 


NACA  TN  2145 


induces  lift  in  region  I*.  The  cancellation  of  lift  in  region  I' 
induces  lift  in  region  I",  and  so  forth.  Each  of  these  cancella¬ 
tions  is  handled  as  previously  described.  These  computations  are 
very  tedious  when  lift  is  induced  upstream  of  a  subsonic  leading 
edge  (for  example,  region  I*  of  fig.  8),  inasmuch  as  a  knowledge 
of  the  sidewash  (Av^.,),  as  well  as  of  the  lift  distribution 

(Aujt),  is  needed  in  order  to  continue  the  cancellation  process. 
Numerical  methods  are  generally  required. 

Successive  cancellations  are  discussed  more  extensively  in 
references  3  and  4. 


SUMMARY  OF  ANALYSIS  AND  APPLICATIONS 

A  general  expression  was  determined  for  the  lift  distribution 
over  a  cancellation  wing.  The  expression  is  valid  when  the  plan- 
form  boundary  (on  cancellation  wing)  separating  the  region  of  zero 
upwash  from  the  region  for  which  the  lift  is  specified  is  everywhere 
subsonically  inclined  to  the  free  stream.  This  expression  permits 
the  determination  of  lift  distributions  on  wings  that  cannot  be 
solved  by  cancellation  techniques  based  on  conical  superposition. 

The  boundary  conditions  for  either  the  flow  about  a  subsonic  leading 
edge  or  a  subsonic  trailing  edge  can  be  satisfied. 

The  lift  cancellation  technique  was  illustrated  for  swept  wings 
having  curvilinear  plan  forms.  Leading-edge  and  trailing-edge 
cancellations  were  considered.  In  addition,  the  loading  in  a  region 
influenced  by  a  reentrant  side  edge  was  found. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  January  16,  1950. 
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APPENDIX  A 

SYMBOLS 

The  following  symbols  are  used  in  this  report: 
(i+e)r-(i-e)s 
(l+aJr^sJ-Cl-Sjs 
(l+0)r2(s)-(l-0)s 

local  lift  coefficient, 
constant 

root  chord  of  swept  wing 


F(0,k) 


Pb-PT 


elliptic  integral  of  second  kind, 
/"'Sin  0 


E(0,k)  = 


l-k^(t) 


d% 


CO 


elliptic  integral  of  first  kind, 
psin  0 


F(0,k)  = 


dttL 


0  ^(i-0)o2)(i-k2«o2) 


G(r,s0)  function  of  r  and  sQ  defined  by  equation  (12) 


2aU 


PE  (§,'fc 


et 


slope  of  plan-form  edge  in  r,s- coordinates,  dr/ds 
modulus  of  elliptic  integrals 
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p  density 

T  area  of  integration 

0  amplitude  of  elliptic  integrals 

<5P  perturbation  velocity  potential 

doublet  strength,  CpT-qpg 

C0o  integration  variable 

Regions : 

I  region  on  cancellation  wing  for  which  loading  is 

specified 

Ia > Ift ,  .  .  .  subdivisions  of  region  I 

II  region  on  cancellation  wing  for  which  w  =  0 
IIa,II^,  .  .  .  subdivisions  of  region  II 

HI  additional  region  on  cancellation  wing  for  which 

w  =  0 

Special  designations: 

r  =  r-^s)  r  as  function  of  s  along  plan-form  boundary  1 

s  =  s^(r)  s  as  function  of  r  along  plan-form  boundary  1 

y  =  yi(x)  y  as  function  of  x  along  plan- form  boundary  1 

x  =  xi(y)  x  as  function  of  y  along  plan- form  boundary  1 

r  =  r2(s)  r  as  function  of  s  along  plan-form  boundary  2 

s  =  s2(r)  s  as  function  of  r  along  plan-form  boundary  2 


and  so  forth. 
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Subscripts: 

1,2,3 

I,II 

B 

T 


refers  to  plan- form  boundaries  1,  2,  and  3,  respectively 
refers  to  regions  I  and  II,  respectively 
bottom  surface  of  z  =  0  plane 
top  surface  of  z  =  0  plane 

variable  of  integration  . 


1297 


NACA  TN  2145 


27 


APPENDIX  B 


» 


SOLUTION  OF  INTEGRAL  EQUATIONS 


Consider  the  following  integral  equation  (in  the  notation  of 
the  appendix  in  reference  5),  where  the  function  f(x)  is  assumed 
known  and  the  function  u(|)  is  to  be  determined: 


f(x) 


ox 


urn  de 


ox 


J 


[uU)-u(x)]  de  2u(x) 

(X-!)3/2  (I^T2 


(Bl) 


After  an  integration  by  parts,  equation  (Bl)  may  be  written 


(B2) 


Equation  (B2)  is  now  an  integral  equation  of  the  Abel  type.  The 
continuous  solution  for  u(|)  is  (reference  14) 


u(z) 


OZ 


1 

2tc 


U  a 


f(x)  dx 


(B3) 


evaluated  at  z  =  (; .  This  result  is  presented  in  reference  5 . 

Equation  (lib)  corresponds  to  equation  (Bl)  with  u(|)  =  G(r,sQ) 

and  f(x)  =  0.  The  solution  for  G(r,s0),  according  to  equa¬ 
tion  (B3),  is  then 


G(r,s0)  =  0 


(13) 
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APPENDIX  C 


DIFFERENTIATION  TO  OBTAIN  Aun 
The  differentiation  indicated  in  equation  (16) 


2J3  .  (  d  ,  d 

“M  AuII  =  V5?  +  Si 


is  to  be  conducted. 
First 


AJ  r-r2(s) 

JC 


pr2(s) 


A<Pj  ^o 


^(s)  (r-roWr2(s)-ro 


(16) 


2p  A 
M  AUH 


1  - 


dr2(s) 


ds 


2ir/^r-r2(s) 


nr2(s) 


ACPl  dro 


>i(s)  (r-ro)A|r2(s)-ro 


^p(s) 


A^Pt  dr 
I  o 


i(g)  (r-r0)/y/r2(s)-r( 
Inasmuch  as  ACp^  is  a  function  of  rQ  and  s. 


n r2(s) 


^o(s) 


S 

c*r 


ACpj  drQ 


ACPj  dr( 


(Cl) 


(r-r0)Ajr2(s)-r( 


r  (s)  (r-ro)Mr2 (s)"r< 


(C2a) 


30  NACA  TS  2145 

and 


so  that  equation  (C2b)  can  be  written  as 
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Upon  substituting  equations  (C2a)  and  (C3)  into  equation  (Cl)  and 
integrating  by  parts  those  integrals  containing  tSP-j-,  equation  (Cl) 
finally  reduces  to 
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APPENDIX  D 


ALTERNATE  DERIVATION  OF  SOLUTION  SATISFYING 


KUTTA  CONDITION  AT  r 


(s) 


The  integral  equation  formulation  in  terms  of  ACp  (equa¬ 
tion  (14))  resulted  in  a  solution  that  was  continuous  in  ACp 
(equation  (15))  but  discontinuous,  in  general,  in  the  derivative 

=  Au  (equation  (17a))  at  r  =  r2(s).  In  order  to  obtain  a 

solution  continuous  in  Au,  an  integral  equation  may  be  formulated 
that  is  similar  to  equation  (14),  but  in  terms  of  Au  rather  than 
ACp.  The  inversion  shown  in  appendix  B  should  result  in  a  solution 
that  is  continuous  in  Au  but  discontinuous  in  the  derivatives  of 
Au  at  r  =  r2(s) . 


Consider  equation  (10)  for  the  w  distribution  in  the 
z  =  0  plane.  This  equation  will  be  differentiated  with  respect  to 
x  using  a  technique  introduced  in  reference  15  (equations  (1)  to 
(3)).  The  expression  for  w  at  any  point  (r,s)  is,  from 
equation  (10), 


ACp  dr  ds 
o  o 


r-r0?72 


(11) 


where  T  is  the  area  abc  in  figure  9(a).  The  wing  is  moved 
upstream  a  distance  <3x  (fig*  9(h) ),  keeping  the  coordinate  system 
fixed  in  space.  The  expression  for  the  upwash  at  (r,s)  now 
becomes 


ds 

o 

(r-r0) 


(D2) 
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The  second  term  on  the  right  side  of  equation  (D2)  is  zero  because 
ACp  =  0  along  the  leading  edge.  Subtraction  of  equation  (Dl)  from 
equation  (D2)  then  yields 


bv  M 

H  “  "  8* 


dso 

(i-r0) 


372 


(D3) 


For  points 
for  the  wing  of 


in  region  II  of  a  cancellation  wing, 
figure  3(a), 


Thus, 


0  = 


dso 


nr 


Au  dr 


rl(so) 


(r-rj 


o 

372 


(M) 


This  equation  is  the  same  as  equation  (lla)  except  that  Au  replaces 
A«P.  The  inversion  by  Abel’s  integral  equation,  for  AuIX  in  terms 

of  Auj  then  gives  (from  equation  (15)) 


Au 


A^r-r  (s) 


II 


r?(s) 


AuT 


rl(s) 


(r-r0)/^2(s)-r0 


(24a) 


Inasmuch  as  the  integral  equations  of  reference  5  are  formulated 
in  terms  of  Au  and  are  inverted  by  means  of  Abel ' s  integral  equa¬ 
tion,  only  solutions  satisfying  the  Kutta  condition  will  result 
therein. 
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(c)  Cancellation  wing 


Figure  1.  -  Superposition  to  obtain  lift  on  given  wing  by 
canceling  lift  on  two-dimensional  wing.  (Given  wing  equals 
two-dimensional  wing  minus  cancellation  wing. ) 


nation  for  obtaining  loading  in  tip  region  of  swept  wing 
having  supersonic  trailing  edges. 
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Figure  6.  -  Cancellation  for  obtaining  loading  in  region  influenced  by  reentrant  side 

edge* 
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Figure  8.  -  Successive  cancellations 


